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Abstract 

Dynamic dielectric and mechanical responses of the ferroclastic-dielectric mi- 
crocomposite were studied in this paper. The mechanical inclusion-matrix in- 
teractions have influence on the mechanical moduli of the composite. We have 
described a mechanical response of the composite which consists of the material 
M and of the other material I with dispersion of mechanic moduli of particles 
due to near sphere shape and due to their anisotropy. Larger dispersion of 
mechanical moduli of the second material leads to softening and the stiffness 
has lower values for small concentrations. The inclusions of the dielectric ma- 
terial in the ferroelastic matrix has the effect also on the dielectric response 
of the microcomposite. We have found under which conditions the inclusions 
have the effect of enhancement of the dielectric response. Results of low fre- 
quency dielectric constant and dielectric loss in orthorhombic ^^2(^^04)3 show 
for polycrystalline material, where voids may play the role of dielectric material, 
linear increasing with increasing hydrostatic pressure. This is qualitatively in 
agreement with our theory. We are explaining the presence of the hysteresis 
phenomena observed in these experiments by amorphization which changes free 
volume and its effective mechanical properties. 



1 Introduction 

The problem of effective sfiear and bulk moduli, of an effective Poison's ratio 
and of an effective dielectric response in microcomposites of the ferroelastic- 
dielectric type is a very interesting one. The static and dynamic dielectric re- 
sponse of ferroelectric-dielectric and the static response of ferroelastic-dielectric 
type microcomposites was studied recently [1] - [6]. A coupling of the elastic 
strain to the electric polarization leads to the dielectric response of a ferroelas- 
tic type material, see in [7] - [8] . A high anisotropy of single crystal properties 
induces an enhancement in the macroscopic piezoelectric response for polycrys- 
talline materials like BaTiOs and PZN - PT, [9]. The authors of the cited 
paper do not consider pores in the material. However there is always a free 
volume in a polycrystalline material which is formed by these pores. We can 
consider polycrystalline materials as a microcomposite in which one type of 
material is a crystalline material and the other one are pores. The authors of 
the mentioned paper determined macroscopic properties of a polycrystal from 
the resulting average strain tensor. They compare results of their modelling 
with experimental data. While they have found that texture and the corre- 
sponding grain size enhance piezoelectric response, one may expect that also 
the anisotropy of clusters increases when the concentration of one component of 
the microcomposite increases leads to similar effects, when the clusters become 
very large. The clusters are not isotropic, they are anisotropic. Thus we may 
expect that in microcomposites of the ferroelastic-dielectric type increase of the 
concentration of dielectric material leads to the enhancement of the dielectric re- 
sponse. Materials of the perovskite type {LaAlO^, CaAlOs, SrAlO^, BaTiO^, 
PbNiOs, Pb{Zr,Ti)03, ...) undergo a phase transition from a cubic phase to 
a phase with lower symmetry at some critical temperature. The dynamic me- 
chanic response, and the dynamic dielectric response of such materials which is 
due to coupling between the elastic strain tensor and the polarization is thus an 
interesting problem. 

We have studied recently the static response of microcomposites of the type 
mentioned above [1], and discussed qualitatively the results of the model for 
pressure dependence of a dielectric constant in the Al2{WO)4 polycrystalline 
material, in which a second material were porous in between crystallites. Let us 
note that the ^4/2(^0)4 type materials are used in many applications like fuel 
cell electrolytes, gas sensors, laser materials, composites with tunable thermal 
expansion, see in [2] and in references therein. Mechanical analysis measure- 
ments are done usually at low frequencies (O.lHz - 10 Hz), however measure- 
ments at higher frequencies are done also, [10] - [11]. Elastic response function 
(compliance) shows Cole-Cole diagrams of the circular and of the non-circular 
behavior for these materials in their crystalline and in their ceramic form [12]. 
Multirelaxation phenomena exist in these materials under some conditions. 

It is known that higher electric and mechanical loading leads to a nonlinear 
behavior in ferroelectric and ferroelastic ceramics [13]. At higher temperatures 



however there is present no ferroelastic phase, there is present a paraelastic 
phase. Couphng of the clastic strain and electric polarization does however 
exist at these temperatures too. Mechanical and dielectric response depend then 
on mechanical forces acting on the microcomposite of the ferroelastic-dielectric 
type at these higher temperatures. Constraints due to neighboring material lead 
to a transition from the paraelastic to a ferroelastic phase [14]. We assume in 
this paper that there are small mechanical fields of the order 10~^ which are 
acting on the microcomposite, and that particles in the microcomposite have 
a diameter of the order of 1/im, see also [1]. For large concentrations of the 
dielectric material the clusters of the ferroelastic phase are small and finite. 
At the percolation transition there are very large clusters of ferroelectric and 
dielectric materials. It is not known how they are responding to the external 
mechanical fields mainly due to the fact that the form of the clusters is irregular. 
To study the region of the percolation transition in the microcomposites of the 
ferroelastic-dielectric type we need to understand the behavior of elastic moduli 
of the microcomposite at this region. 

In a ferroelastic phase long-range anisotropic forces may appear [15] in the 
region when large inclusions are present. In our paper we will discuss proper- 
ties of the microcomposite above the critical temperature. Thus the long-range 
forces can be neglected. To determine the behavior of elastic moduli of the 
microcomposite we assume that the grains of both phases are more-less of the 
spherical shape. The shape of the ferroelastic material particles is assumed 
to be spherical and characterized by the same diameter, it dependens on the 
method of preparation of the microcomposite. The shape of dielectric particles 
is assumed to be more-less spherical and we assume further that there exists a 
distribution of their radii and of their shapes which may be characterized by a 
homogeneous distribution function as concerning mechanical properties. Analy- 
sis of the elastic moduli of heterogeneous materials was done in [16]. We will use 
this method to analyze our system of the ferroelastic-dielectric microcomposite. 
We will find mean-field equations for the effective shear and bulk moduli of the 
composite, and for the effective Poisson's ratio. 

The aim of this paper is to study the dielectric and mechanical response of 
microcomposites of the ferrolastic - dielectric type. The mean-field equations 
for the effective shear and bulk moduli of the microcomposite together with 
the mean field equation for the dielectric susceptibility enable us to study these 
responses. 

2 Model of ferroelastic-dielectric microcompos- 
ites 

Let us describe a model for ferroelastic-dielectric microcomposites. The po- 
larization, a primary order parameter in fcrroelectrics, is in ferroelastics a sec- 



ondary order parameter. Primary order parameter in these later materials is the 
corresponding component of the elastic strain tensor (or a combination of com- 
ponents of the elastic strain tensor) . The ferroelectric state is present in such a 
material due to a coupling between the elastic strain tensor and the polarization. 
Changing the concentration of two types of particles in the microcomposite the 
response to external fields changes. This holds for dielectric response and for 
mechanical responses. As we noted for simplicity we consider all particles in 
the microcomposite of the same diameter d for the ferroelastic material and of 
nearly the same diameter and of nearly spherical shape for dielectric material, 
which is a realistic approximation. This enables us to study microcomposites in 
which dielectric particles may be also porous, and thus polycrystaline materials. 
The concentration at which an infinite cluster of dielectric particles forms is the 
percolation transition concentration. Increasing further the concentration of 
the dielectric material this material becomes the matrix material. Thus the in- 
terplay between the percolation transition and the ferroelastic phase transition 
appears. 

3 Elastic moduli of the microcomposite 

We will assume that the microcomposite consists of particles - spheres of the 
material M, we will use the index N for the material constants of this material, 
and of particles - almost spheres which have different " diameters" of the other 
material. We will assume further that there are (N-1) particles of the second 
type. Mechanically we thus have a composite from N types of particles. Let 
us assume that the concentration of the particles of the M material is 1-x, 
and that the i = I, ...,N — 1 types of particles of the second material I have 
the same "concentration" Xi = jf^ . Thus the concentration of the second 
material is x. We can imagine this composite as if it is to consisted of regular 
grains of the material M and of less regular, even almost irregular, grains of 
the second material which are almost spheres and which have some distribution 
of diameters. The size distribution of grains of the second material is assumed 
to be homogeneous. Then mechanically we have a material of the volume V 
with the N phases. The volume concentration is -^ for the i-th particle. We 
assume that the particles of the second material have the same volume -^^ . 
Derivation of effective moduli using the method according to [16] leads to the 
effective shear modulus G* and to the effective bulk modulus K* and to the 
effective Poissons ratio ly* . We consider, as in [16], a large cube with edges 
parallel to the coordinate axes (x,y,z). We further assume that a uniform shear 
stress acting to the surface of the cube and a uniform hydrostatic pressure acting 
on the surface of the cube lead to a strain and a volume contraction. Note that 
for N large the particles of the second material in the material M will appear 
as inclusions in a matrix consisting of that N phase. However in [16] it is used 
the composite as a matrix, and thus he includes interactions of these particles 



of the second material in the material M as the concentration of these particles 
increases. This then leads to effective elastic moduli. 

Let us now describe equations for the effective elastic moduli in the limit of 
a large number of particles of the second material. The effective shear modulus 
G* is found from the equation: 

G*- Gn^'^Jo ^^ GN'G*+l3*{G{r)-G*) ^^' 

where: 

'^ 15(1-1/*) ^' 

and v* is the effective Poissons ratio of the composite material. 
The effective bulk modulus K* is found from the equation: 

1 _ 1 f^ K{r) 1 



where: 
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3(1 -z/*) 
and V* , the effective Poissons ratio of the composite material, is given by: 

^ ,K*-2G* 

QK* + 2G* ^ ' 

Thus equations (1) - (5) give the dependence of G* and K* on the elastic 
constants Gat and K^^^ and G{r) and K{r), the distribution functions of shear 
and bulk moduli for the second material. Here r denotes different types of 
particles of the second material. 

Let us assume that the distribution of the clastic moduli for particles of the 
second material is linear: 

G{r) = rGi + (1 - r)Go (6) 

and: 

K{r) = rKi + (1 - r)K„ (7) 

where Go and Gi are two values of the shear modulus G for the second 
material which correspond to the index and 1 of the particles, and where 
Kq and Ki are two values of the bulk modulus K for the second material which 
correspond to the index and 1 of the particles. Let us note that the distribution 
may be in general different from the linear, however the linear approximation is 
expected to be good for the case in which values for the moduli for r = and for 



r = 1 are not different too much. Note that for larger difference Gi — Go , i.e. 
for the second material in which rigidity of particles changes with the change of 
the shape and the size, e.g. due to presence of texture, cracks or other types of 
defects. Then we obtain equations (1) - (5)in the following form. 
The effective shear modulus G* is found from the equation: 



±^J_ X P*{G*~Gn)-G* G*+(3*{G,-G*) _ 

G* Gn Gn(3* /3*(Go-Gi) '^^G*+f3*{Go-G*)' ^ ^' 

where /3* is given (2) and u* is again the Poissons ratio of the composite 
material. 

The effective bulk modulus K* is found from the equation: 



K* Kn KNa*^ a*{Ko-Ki) ^ K* + a*{Ko - K*)^ ^ ^' 

where a* is given in (4). 

From the equations (8) - (9)it follows that such textures, cracks and defects 
which induce larger difference of the elastic moduli G and K lead in the first 
approximation of the inverse of the absolute value of the difference Gi — Gq and 
of the difference Ki — Kq to the microcomposite in which: 

1 1 . a; , , ^ 

and: 

and where for a* and /3* positive hardening of the microcomposite is present 
in this order. This is observed in numerical simulation, [9], where however also 
crystallographic microstructure is taken into account in the numerical modeling. 

The equations (8) - (9) have different form for the case in which there is no 
distribution of the clastic moduli of particles of the second material, e.i. when 
Go = Gi = G and when i^o = Ki = K . Then we obtain that the equations (8) 
- (9) take the following form. 

The effective shear modulus G* is found from the equation: 

± = J ^ r ^-^^ 1 (121 

G* Gn Gn G*+f3*{G-G*)' ^ ' 

and the effective bulk modulus K* is found from the equation: 
1 1 X ^ K-Kn 



K* Kn Kn K* +a*{K -K* 



(13) 



In this chapter we have formulated a theory for the elastic moduli of a micro- 
composite which has two constituents, each of them is isotropic and elastic, the 
first one has the diameter of particles the same, the second one has a distribution 
of the form of particles which are almost spheres and which have distribution of 
their elastic moduli, its concentration is x. We studied in (12) - (13) the case 
of no distribution of the form of particles of the second type. 

4 The case of no distribution of the form of par- 
ticles and effective elastic moduli 

In (12) - (13) we have described the equations for the effective shear and bulk 
moduli in the case of no distribution of the form of particles of the second 
type. To study the solution of these equations it is convenient to introduce the 
following ansatz in the equations (8) and (9): 

G* = {G-Gn).9 (14) 

K* = {K-KN).k 
We have found that the functions g and k arc given by: 

where: 

Dg = r (7 + in) -{x + iN)f + 4(1 - nP*llN (16) 



and where: 



and: 



IN = -^ — -^ (17) 
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and: 



where: 
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k± = ^^ j.ha*(<5 + <5w) + (,x + 5Ar)±v^D^] (19) 

z[l — a \ 



Dk = [a*{S + 5n) -{x + Sn)^ + 4(1 - a*)a*SSN (20) 



and where: 



and: 



^-X^ (-) 



-^. (22) 

K-Kn ^ ' 



The effective Poison's ratio i'* has the form: 

^ , ^ 3fc - 2ge 
6fc + 25e 
where e is defined as: 



(23) 



Note that Dg and I?^ are discriminants which may be positive or negative. 
In the first case g± is real and k± is real too. In the second case g± is a number 
and k± is a complex number too. Note also that the signs of both discriminants 
in equations (27) and (17) are independent. Different signs may be valid in 
different regions of concentration. Here we will consider the static case only in 
the following. According to the fluctuation-dissipation theorem the imaginary 
part of the moduli will be zero. Neverthless note that the results above are 
obtained for general moduli and that complex moduli can be measured. 

For isotropic elastic solids [17] the range of i^* of stability is between -1 and 
0.5. We will assume that the microcomposite is an isotropic elastic solid. Thus 
if a percolation transition in this microcomposite occurs at ly* > 0.5 it is behind 
the stability region. Then increasing the concentration of the second compo- 
nent /or decreasing concentration of the second component if this component 
forms a matrix/ may lead to an unstable microcomposite. The microcomposite 
is assumed to be unconstrained. Such a results would predict, if we assume 
that the isotropic elastic solid /microcomposite/ remains isotropic changing the 
concentration, that it may become mechanically unstable. However let us note 
that in [17] the composite as a continuum isotropic is considered, and this may 
be not true for large clusters. 

The range of the Poison's ration v* above corresponds to the stability crite- 
ria, e.i. to positive G* and K* mechanical moduli. Note however that ferroelas- 
tic materials may have a negative stiffness, see in [18]. Then the problem with 
negative stiffness inclusions /clusters/ in a matrix may be considered by our 
model, this problem is an interesting one but not discussed here. Note that in 
[19] ferroelastic martensites parent to product phase transition was studied. Or- 
der parameter which drives this transition is shear deformation. Strain energy 
dominates morphology, and the transition results in a characteristic lamellar or 
twinned structure. This then may lead to a nonisotropic microcomposite, and 



calculations of the mechanical moduli should take this fact into account. It is 
not the case in this paper for simplicity. Martensites are proper ferroelastics, 
and improper ferroelastics are f.e. ferroelectrics, the martensites show weakly 
first order phase transition near the second order, or second order phase transi- 
tion and thus Landau theory applies. This fact is used in our paper. Examples 
of proper ferroelastics are InTl, FePd, NiTi, AuCd. In proper ferroelastics 
textures form occurs and it does not correspond with the assumption of the 
isotropic elastic microcomposite. Improper ferroelastics are materials like high- 
temperature superconductors, and GMR manganitcs. In [20] it was studied how 
composites with inclusions of negative bulk modulus in a matrix behave. Let 
us note that complex values of mechanical moduli correspond to viscoclastic 
materials, in which however G , the imaginary part of the shear modulus, must 
be positive [21] if the material should be stable. The quantity ta,n{5) is given 

by the ratio ^, where G is the real part of the shear modulus. Expcrimen- 
tal results for this ratio arc from 0.01 at room temperature to 0.015 at higher 
temperatures. In [22] it was found that riples arc present in the response for 
temperature dependence of the 5-percent VO2 in Sn matrix. 

The shear stress induces large changes in strain when the shear modulus 
is small in its value. At the percolation transition when one cluster of the 
dielectric material becomes infinitely large cluster, the mechanical properties 
of the microcomposite may be expected to be such that small stress induces 
large changes in the strain, even that at this transition and near this transition 
an instability occurs. The material becomes very susceptible to the stress, or 
even mechanically unstable. While pressure induces changes in clusters of both 
types of materials in the microcomposite, and thus the susceptibility of the 
microcomposite to pressure will be present too, the clusters change their volume 
while interacting, the shear stress leads to shear stress of clusters acting on other 
clusters. The clusters do not have a regular shape. Due to this fact we may 
expect that pressure on the microcomposite leads to shear stresses between 
clusters inside the material. And the shear stress acting on the microcomposite 
leads to volume changes of clusters induced by other clusters, on which the shear 
stress is acting. Thus at and near the percolation transition we may expect that 
both mechanical moduli become /almost/ zero. 

We have found that expansion in (3* is useful. When (3* tends to zero then 
G* tends to small values. This is found to occur for v* — ^, a* is finite. The 
microcomposite is however unstable under these conditions. Zero (3* would cor- 
respond to a critical concentration Xc which may be calculated. This is not so 
unnatural because the clusters change their size going with the concentration to 
the percolation threshold, thus volume properties of the components are chang- 
ing depending on their mechanical properties. Note however that the v* = ^ 
value is beyond the stability region of the material, as noted above. So there 
exists another critical value of the concentration, a;„i, above which the compos- 
ites become mechanically unstable depending on their material constants. We 
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described these results qualitatively, it is not the aim of this paper to study 
mechanical instability of microcomposites which may appear for certain range 
of material constants. 



5 Low concentration of the dielectric material: 
elastic moduli 

From the equation (12) - (13) we obtain in the limit of small concentrations x 
up to 0{x) of the dielectric particles the effective Poison's ratio: 

, 'SKn - 2Gn , „. . . 

where the constant is given by: 

^ ^ iKN - 2Gn 6Kn + 2Gn 

and the effective bulk modulus: 

K*^K^{1--^—) (27) 



Sn + a*S' 



and the effective shear modulus: 



G*^GNil- ) (28) 

In + P 7 

where a and [3 constants are calculated as above for the material M. 
We have found that: 

and for a* we have found: 

a =aAr(l-x(— — T + iTf. r) (30) 

i{l + VN) ?>{1-vn) 

Note however that in the equations (27) - (28) we can omit the star at the 
quantities a and fi due to the fact that the calculations are done up to the 
second order in the concentration x, i.e. up to the quantities of the order 0{x'^). 

From the equations (27) we can find that the material with inclusions is less 
harder if the microcomposites is such that the quantity 5n + aw (5 is positive. 
And vice versa. From the equation (28) we can find that the material with 
inclusions stiffens if the microcomposite is such that the inequality ^m + aN"t 
is positive. And vice versa. 
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6 Low concentration of the dielectric material: 
dielectric response 

For low concentration of dielectric particles in the ferroelastic matrix we may 
calculate the dielectric response of the microcomposite as in [1] . For values of 
small X the effective dielectric permitivity e^ff is given by: 

aj afLrM 

1 QTpGn , , 1 ,, , QTp , ^d-ij{l + ^^) 
I o, I 2,x — (1 I > J J 



a 



Sm+umS a/' afGN \d + 2-^{l + -^^ 



where a/ is the dielectric permitivity of the ferroelastic matrix. We see 
from the equation (31) that when the material elastically softens due to increas- 
ing concentration of dielectric material, i.e. when the conditions for softening 
described above are fulfilled, then the response of the ferroelastic-dielectric mi- 
crocomposites is higher. This is true for small concentrations of the dielectric 
material. 

7 Summary 

We studied here dielectric and mechanical responses of the microcompositcs: 
ferroelastic-dielectric. A model for such a microcomposite was formulated in 
[1]. Cubic symmetry of the material was assumed which corresponds to most of 
the known ferroelastic materials at high temperatures. Nevertheless our results 
give a qualitative picture also for other symmetries. In our paper [1] we did not 
consider the mechanical inclusion-matrix interactions. This has influence on 
the mechanical moduli of the composite. The ferroelastic crystals as ferroelastic 
particles are assumed here to have no domain wall structure. We have described 
a mechanical response of the composite which consists of the material M and of 
the other material with dispersion of mechanic moduli of particles due to near 
sphere shape and due to anisotropy. Larger dispersion of mechanical moduli 
of the second material leads to softening and the stiffness has lower values for 
small concentrations. For no dispersion case we have found a known result 
for the effective moduli of the microcomposite. The inclusions of the dielectric 
material in the ferroelastic matrix have the effect on the dielectric response 
of the microcomposite. We have found under which conditions the inclusions 
enhance the dielectric response. Effect of elastic clamping was not discussed 
in this paper. For its discussion for improper and pseudoproper ferroelastic 
inclusions see [23]. 

Results of low frequency dielectric constant and dielectric loss on orthorhom- 
bic ^4/2(^1^04)3 show for polycrystalline material [2], where voids may play the 
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role of dielectric material, linear increasing dependence with hydrostatic pres- 
sure. This is qualitatively in agreement with our theory. The authors of these 
measurements observed hysteresis with increasing and decreasing pressure. This 
hysteresis may be qualitatively understood within our model due to amorphiza- 
tion of the polycrystalline material at higher pressures and due to the fact that 
the aniorphization remains decreasing then the pressure. Amorphous micro- 
composite has different mechanical properties than polycrystalline. From our 
theory it follows that for different mechanical constant values there are different 
slopes and dielectric response dependence on the pressure. In epitaxial films 
proper ferroelastic phase transitions with symmetry-conserving and symmetry- 
breaking misfit strains may be present [24] . Most of microcomposites arc in the 
form of films. The authors [24] have found that if the extrinsic misfit strain 
does not break the symmetry of the high-temperature phase, the transition in 
the film occurs at somewhat lower temperature than in the bulk. This may play 
the role in the observation of the transition in this material. 

In cubic-tetragonal systems like Nb^Sn, V^Si, In — Tl alloys, Fe — Pd alloys 
and Ni2MnGa the cubic cell elongates (or contracts) in one of the main axis 
to form tetragonal cell [25]. While these interesting materials are in general 
different from those about which we wrote in the Introduction and above, their 
high-temperature phase is cubic and our approach may be used to study them. 
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